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ABSTRACT 

The thermal disturbance observed during the sudden stratospheric warming 
in January 1963 showed a wave -like behavior. The wave characteristics of this 
event were consistent with the propagation and attenuation of downward propa- 
gating turbulent heat conduction waves . 

Full wave calculations indicate that the stratosphere behaves like a selec- 
tive filter with respect to atmospheric waves. Waves can propagate northward 
only under the condition of a resultant southerly wind of about 5 m sec" 1 . Up- 
ward propagating waves in the mesosphere are attenuated in such a manner that 
their relative amplitude generally does not increase with altitude. These con- 
siderations, in conjunction with the physical interpretation of observations of 
sudden stratospheric warmings and simultaneous short-period winter anomalies 
in the ionospheric D-region, lead to the suggestion of a common source, in form 
of a sudden initial heat impulse, of both phenomena. An upper limit for the 
generation of the initial heat impulse was placed at 50 km on the basis of ob- 
served wave characteristics and of static stability requirements in the polar 
winter mesosphere. Thus, the common source seems to be located at or a little 
below the stratopause. 


*NAS-NRC Research Associate on leave from the Astronomical Institutes of the University of 
Bonn, Germany 

f NAS-NRC Research Associate on leave from the Institute for Meteorology and Geophysics, Free 
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' YRATOSPHERIO-MESOSPHERIC COUPLING 
DURING STRATOSPHERIC WARMING 

1. INTRODUCTION 

From investigations by BOSSOLASCO and ELENA (1963) and others (SHAPLEY 
and BEYNON, 1965; LAUTER, 1967; BELROSE, 1967) it is now evident that sudden 
stratospheric warmings during winter, discovered by SCHERHAG (1952) show 
significant correlation with the increase of nondeviative absorption of medium 
and high frequency radio waves in the ionospheric D-region, known as winter 
anomaly of absorption, and with the simultaneous decrease of the virtual reflec- 
tion height of LF- and VLF-waves. These ionospheric effects are believed to be 
closely related to changes in the neutral air properties, like temperature and 
pressure, and (or) variations of photochemical processes due to changes in 
large-scale neutral wind circulation during breakdown of the stratospheric 
polar vortex in the winter hemisphere (SECHRIST, 1967), 

There remains the question of the origin of this disturbance and of the trans- 
port mechanism within the atmosphere. REGENER (1952) suggested already that 
the trigger mechanism for the stratospheric warming phenomenon may be located 
within the stratosphere, i.e., well below the ionospheric D-region, and FORTAK 
(1958) proposed a vertical turbulent heat transport as the propagation mechanism 
of the disturbance. These two suggestions will be considered in some detail in 
this paper. The wave structure of the sudden stratospheric warming phenomenon 
will be discussed quantitatively in terms of the event of January 1963. It 
will be shown that the disturbance originates from the region of higher tempera- 
ture near the low latitude stratopause or even from the summer hemisphere, the 



easterly circulation of which seems to expand abruptly toward the winter pole as 
a result of instabilities of the general wind system. Full wave calculations of 
turbulent atmospheric waves reveal that the characteristics of heat conduction 
waves are consistent with the observed wave structure during a sudden strato- 
spheric warming. 

2. TEMPORAL BEHAVIOR OF THE 1963 STRATOSPHERIC WARMING 

One of the strongest stratospheric warmings ever observed occurred in 
January 1963. The synoptic features of this event have been described by FINGER 
and TEWELES (1964). The associated energy processes have been analyzed by 
JULIAN and LABITZKE (1965), PERRY (1967), and HIROTA (1967). This 1963 
event was typical for strong stratospheric warmings. Temperature cross- 
sections and time-sections appeared to be very similar to the warmings in 
February 1952 (SCHERHAG, 1952) and January 1968 (FINGER and McINTURFF, 
1968). 

Figures la and lb show temperature cross-sections along 90°W and 90°E 
during the maximum of stratospheric warming effect, on 26 January 1963, and 
six days earlier, on 20 January 1963, respectively. The data were taken and 
analyzed from a series of daily stratospheric maps published by the Institute 
for Meteorology and Geophysics, Free University of Berlin (METEOROLOGISCHE 
ABHANDLUNGEN, 1963-1966) and from WARNECKE (1968). Although the iso- 
thermal planes of the disturbance are tilted not only toward the north but also 
toward the west and propagate toward the northwest, changing the direction of 
propagation with time, the tilt and the direction of propagation are predominantly 
northward. Therefore, we consider a meridional cross section only. Figure 2 
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shows the stratospheric temperature variation with time for a number of isobaric 
levels above the location 70°N, 90°W, near the center of maximum warming. 

From this figure, the time of occurrence of the maximum temperature as a func- 
tion of isobaric level was determined and plotted, Figure 3, versus height (filled- 
in circles). The horizontal solid bars indicate the estimated range of uncertainty 
of the observations. The open circles and the horizontal dashed bars have been 
derived from the observations of the Berlin Warming in February 1952 (FOHTA.K, 
1958). 

In Figure 4 the decrease of the relative amplitude AT/T, normalized by the 
lOmb value, has been plotted versus height. Here T (in degrees Kelvin) is a 
reference temperature, defined as the average of the observed equilibrium tem- 
peratures before and after the warming phenomenon. AT is the maximum devia- 
tion from the reference temperature T. Since the equilibrium temperatures 
before and after the warming phenomenon are not very well defined in the obser- 
vations, the range of estimated uncertainty of the derived parameters is indicated 
by solid bars in Figure 4. The open circles and the dashed bars again represent 
data derived from the 1952 warming. Since this event has been observed over 
Berlin, Germany, i.e., at one random location, only, while the data of the 1963 
warming were taken at the well defined location of maximum effect, some pre- 
caution is advisable in comparing the absolute values of the parameters derived 
for the two events. But nevertheless, very similar features can be seen in 
Figures 3 and 4. The larger time delay in the occurrence of the tempera- 
ture maximum as well as a stronger decrease of the temperature amplitude with 
decreasing height was observed in case of the 1952 warming. 
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3. WAVE PARAMETERS DERIVED FROM TEMPERATURE OBSERVATIONS 

The history of the 19(33 stratospheric warming at different altitudes, as 
demonstrated in Figure 1, suggests that this phenomenon be treated theoretically 
as wave-like propagation of an impulse type disturbance. The main task pursued 
in this study was to determine the wave mode involved from the observations. 

As a first approach, the discussion was confined to the highly idealistic model 
of plane free internal harmonic waves extended infinitely in the horizontal and 
vertical directions. Such waves can be described by the expression 


t = e o e 


jcot - jk x x - j k z z + fa 


where t is time, w is angular frequency. 


CO 


(1) 


(2) 


is the horizontal wave number which is assumed to be real and constant, V is the 

V X 

horizontal phase velocity. A real and constant value of k x is suggested from the 
temporal and spatial behavior of the observed disturbance, which propagates with 
fairly constant horizontal velocity and with only slight changes of its shape (see 
Figures la and lb). 


M z ) = k zr ~ i k zi ( 3 > 

is a complex vertical wave number, with j = /^r. 

k ( Z ) ~ (4) 

zr ^ ' V z ( z ) K 
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is the real part of k z . The altitude dependence of k 2r involves diffraction of the 
waves and therefore a bending of the wave planes. V z is the vertical phase velocity. 
k zL (z) is the negative imaginary part of k x , responsible for the dissipation of 
wave energy in the direction of propagation. Therefore k zi is positive for waves 
propagating upwards (transporting energy upwards) and is negative for waves 
propagating downward. The wave amplitudes of our wave modes ax*e normalized 
in such a manner that their squared magnitude is a measure of wave energy (see 
Appendix B). Therefore a term 1/2H has been added in the exponent in Equa- 
tion (1), where H is the scale height. This takes into account the adiabatic 
increase with height of the amplitudes of physical parameters like relative 
temperature disturbance AT/T or relative pressure disturbance Ap/p within an 
isothermal atmosphere. Thus 


1 

2H 


<1 f u Pgoing waves 

k .l > 2H for • 

dovvngomg waves 


(5) 


may be larger or smaller than zero for upgoing waves indicating an increase or 
decrease of the relative amplitudes with height. Downward propagating waves 
always show a decrease of amplitude in the direction of wave propagation. 


According to Figure 1, the actual wave structure is considerably more com- 
plicated than assumed in the plane wave model. In fact, the horizontal and 
vertical velocities of the disturbance determined from Figure 1 are group 
velocities rather than phase velocities. Furthermore, the wave is attenuated 
not only in the vertical but also in the horizontal direction, suggesting spherical 
waves rather than plane waves. Therefore, only qualitative agreement between 
the wave characteristics derived from Figure 1 and the theoretical model, 
outlined in section 4, can be expected. 
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From Figures 1 to 4 the following wave characteristics were determined: 

(a) components of the mean group velocities, supposed to be nearly equal to the 
phase velocities 


V x - 4,5 x 10’ 3 km sec" 1 

( 6 ) 

V z ~ 2.5 x 10" 5 km sec" 1 ; 

(b) mean angle of incidence: 

V z 

0 = arctan-rjr*~ Of 32; (7) 

v x 

(c) the predominant angular frequency of the impulse: 

2rr .. _ , 

w ~ ~ 7 x 10"° sec 1 (8) 

where At ~ 10 days is the width of the impulse. 

From Equations (2), (4), (6), and (8) it follows that 

k x ^ 1.5 x 10" 3 km" 1 

(9) 

k ~ 0.28 km" 1 

z r 

For comparison, the westward tilt given by HIHOTA T s (1967b) Figures 3c and 
3d along the 50°N latitude circle is 


0 - Of 13 
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and the horizontal west component of the velocity of the disturbance is 

V x ~ 2. 1 x 10" 3 km/see , 

These data are smaller than the values of the north components by a factor of 
more than 2 and become even smaller near the center of the disturbance. This 
justifies neglecting the western component of the disturbance in our calcula- 
tions. 

From Figure 4 and Equation (5) we find the mean value of the attenuation 
rate as 


~ -0.05 knT 1 . (10) 

The dashed line in Figure 4 is related to the amplification of waves according to 
the term l/2 H. Since the observed attenuation rate is negative, and therefore 



( 11 ) 


it can immediately be concluded, from Equations (5) and (11), that the wave 
propagates downward, transporting energy in the downward direction. 

4. ORIGIN OF THE HEAT IMPULSE 

A complete theory of the sudden stratospheric warming phenomenon has to 
explain not only the wave structure, discussed in the previous sections, but also 
the origin, shape, and direction of propagation of the heat impulse. Due to the 
lack of comparative observational data, the quantitative investigation in this paper 
was mainly restricted to propagation characteristics of the predominant wave. 



With regard to the spatial origin of the heat impulse, only a more or less qualita- 
tive estimate on the initial altitude was possible on the basis of the discussed 
wave theory and sparse observations. 

SHAPLEY and BEYNQN (1965) found no time lag between stratospheric 
temperature changes at 31 km and simultaneous D-region absorption events. 

They suggested a common source for both phenomena at an intermediate altitude. 
LAIJTER (1967) however, found that the D-region abnormal absorption preceeds 
the maximum of the stratospheric warming event by one or two days. From 
analysis of significantly correlated stratospheric and ionospheric data obtained 
over Canada, MANSON (1968) concluded that the origin of the disturbance is 
above 25 Ion. The data presented in Figure 1 of this paper and the interpretation 
of the warming as a wave-like phenomenon suggest an origin at or above the top 
level of the presented observational data, i.e., at or above the 10 mb level or 
31 km altitude. In the following, an estimate of the height of origin of the initial 
heat impulse on physical grounds will be attempted. 


Considering a downward propagating attenuated wave and a lower limit of 
the decrease of relative amplitude of temperature with decreasing altitude, it is 


AT i > AT | -J Z Z Q 
T zo T z 


( 12 ) 


according to Equation (5). Since during stratospheric warmings temperature in- 
creases of 50°K have been observed near 30 km (WARNECKE, 1962), the vertical 
distribution of the lower limit of relative temperature amplitude, represented in 
curve A of Figure 5, was based on this value for 30 km altitude. This value will 
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certainly be on the low side because of the incompleteness of observational data. 

It can be deduced from this figure that a minimal initial temperature impulse of 
100°K at 40 km, 185°K at 50 km, or 370°K at 60 km were necessary to explain the 
observed amplitude of at least 50°I< at 30 km. These altitudes are, however, 
beyond the top altitude of present routine networks of sufficiently dense obser- 
vations in space and time. Thus no complete survey on actual absolute tempera- 
ture amplitudes above 35 kilometers, associated with sudden stratospheric 
warmings, is available for an exact determination of the vertical and horizontal 
localization of the initial heut impulse. A temperature increase of nearly 40°K 
was observed by SECHRIST et al. (1968) between 40 km and 50 km and between 
75 and 85 km at Wallops Island, Virginia, during rapid changes in D-region 
absorption. JONES et al. (1959) reported a temperature increase of 58°K within 
four days observed at 40 km over Churchill, Manitoba. During the recent 1967/68 
stratospheric warming (JOHNSON, 1968) British rocket measurements over West 
Geirinish, Hebrides revealed a maximum temperature change of nearly 90°K 
around 43 km altitude from the warming peak toward the post-warming tempera- 
ture level (LABITZKE and SCHWENTEK, 1968). Unfortunately, no measurements 
were taken during the pre-warming period. 

These observed values, even the latter one, may, however, still be far from 
the total maximum temperature amplitude during sudden stratospheric warmings, 
because in all three cases neither the entire warming period nor the region of 
maximum warming was covered by the rather sporadic observations. Thus, 
these data should be considered more as minimal values. 

An upper physical limit for a temperature increase in the lower mesosphere 
or at the stratopause level is given by hydrostatic stability conditions, i.e., by 
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the dry adiabatic lapse rate which cannot be exceeded. Assuming an isothermal 
mesosphere, which is fairly realistic during winter over higher mid-latitudes 
between 50 km and 70 km (SMITH, et al., 1967), a temperature increase by 200°K 
at 50 km altitude is permissible before vertical instability throughout the major 
parts of the mesosphere occurs. This is illustrated by curve B in Figure 5 which 
represents the permissible temperature increase as a function of height in 
an atmosphere of zero vertical lapse rate between 50 and 70 km. Below 50 km 
allowance was made in curve B for a subsidence of the stratopause during the 
warming process, i.e., the permissible temperature rise increases more strongly 
below the stratopause because of the lower temperatures in the stratosphere. 

As curve B represents maximum values , curve A minimum values , one 
basic conclusion can be drawn from Figure 5: the altitude of origin of the initial 
heat impulse cannot be higher than 50 kilometers, the altitude of the cross-over 
point of both curves, and the height of the origin as well as the initial amplitude of 
the heat impulse can very likely be expected within the hatched region of Figure 5. 

While the vertical localization, of the origin of the heat impulse could at least 
be estimated by physical reasoning, a discussion of its geographic origin and of 
its physical nature would become rather speculative. 

From a meridional temperature cross-section it can be learned that the 
layer of high temperature, centered at the stratopause and caused by absorption 
of solar ultraviolet radiation by atmospheric ozone, fades out toward the winter 
pole and does not exist during polar night poleward of 65° latitude (COLE and 
CANTOR, 1963). During a sudden stratospheric warming the polar vortex, which 
prevents a direct (advective) latitudinal stratospheric heat exchange within the 
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stratopause layer, breaks down. This break-down presently is assumed to be 
caused by baroclinic instability processes within an irregularly shaped vortex 
(MATSUN'O and HIROTA, 1966). But barotropic processes were also suggested 
(JULIAN and LABITZKE, 1965) to accomplish the accumulation of energy available 
for sudden release into sensible heat during the warming event. Recent numerical 
experiments with a rather sophisticated multi-layer atmospheric model by 
MIYAKODA et al. (1968) were successful in simulating the split of the polar vortex, 
which usually is associated with sudden warmings, but failed to produce the sud- 
den stratospheric warming itself out of the baroclinic model. The MIYAKODA 
experiment also showed that barotropic processes were insignificant for the 
polar vortex break-down and the warming phenomenon as well. 

Whether the failure of the simulation experiment was due to residual incom- 
pleteness of the atmospheric model used or the sudden stratospheric warming 
phenomenon is predominantly caused by a yet undetected process or group of 
processes could not be solved. The MIYAKODA Experiment showed, however, 
that barotropic processes were insignificant for the polar vortex break-down and 
the warming phenomenon as well. Thus the process of generation of the heat 
impulse is still unknown, HIROTA, 1967b pointed out that the downward propaga- 
tion of the disturbance during the 1963 warming was preceded by an upward 
motion of a warming process from the 100 mb or 50 mb level into the 10 mb 
region. He, therefore, distinguishes between an intensification stage, during which 
the disturbance propagates upward from the lower into the upper stratosphere, 
and a migration stage, during which the disturbance propagates toward the pole 
and downward, in connection with the break-down of the polar vortex. The 
first stage, between 11 and 15 January 1963 in his case, might act as a 
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trigger mechanism for the second stage, but is certainly not a thermal source of 
the heat energy observed during the migration stage because its heat energy con- 
tent is only a fraction of the heat energy exhibited during the second stage [see 
for example our Figure 2 or IIIROTA’s (1967b) Figure 5] . 

Because of these uncertainties in the present understanding of the generation 
mechanism of the initial heat impulse during sudden stratospheric warmings, this 
investigation was concentrated on the second stage, in HIROTA’s sense, i.e., the 
migration stage, which so vehemently performs transport of heat energy toward 
the winter pole and represents one of the most drastic mechanisms of heat ex- 
change between the sunlit and the polar night sections of the earth’s atmosphere. 
This energy flow has an impulse shape due to the finite abundance of thermal 
energy within the stratopause region. We consider this energy as the generator 
of a spectrum of atmospheric wave modes which propagate in all possible direc- 
tions. The response of the stratosphere and the mesosphere to these waves, 
depending on frequency and angle of incidence, will be discussed in the following 
sections. 

5. MECHANISM OF WAVE PROPAGATION WITHIN THE STRATOSPHERE 

The atmosphere below 100 km is in a highly turbulent state. This can be 
described in a first order approximation by a frequency dependent turbulent 
diffusion coefficient. For periods of more than a day, the turbulent diffusion 
coefficient within the stratosphere and mesosphere is about 10 3 times larger 
than the molecular diffusion coefficient of air (KELLOGG, 1964; NEWELL, 

1963). Eddy viscosity and eddy heat conduction are linearly related to the eddy 
diffusion coefficient. If we assume the factors of proportionality between these 
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three coefficients to be of the same order of magnitude as the corresponding 
molecular coefficients, the coefficient of eddy heat conductivity within the strato- 
sphere becomes 

K eddy ~ A |/tF~ 3.5 x 10 6 erg cnT 1 sec"* 1 deg” 1 (13) 


with A = 2 x 10 s erg cm” 1 sec” 1 deg~ 3/2 and an eddy coefficient of viscosity 
of 


eddy 


8 x 10 


-*/F = 


0. 14 g cm” 1 sec 


1 


(14) 


The exact numerical value of k eddy is, however, not very critical for the following 
discussion, as will be shown later. 

It was shown by VOLLAND (1968a) that in the treatment of the propagation of 
plane internal harmonic waves the coefficient of viscosity can be neglected as 
compared with the coefficient of heat conductivity at very low frequencies [see 
for example, the condition in Equation (A-l)] . This leads to the predominance 
of two wave types in the atmosphere: gravity waves and heat conduction waves. 
The dispersion formula of these two waves can be found in Appendix A, Equa- 
tion (A-2). Both wave types are coupled with each other, in a realistic atmos- 
phere. In order to find the wave pattern of these waves and their relative 
importance, a full wave calculation must be performed. Some basic information 
about full wave analysis is outlined in Appendix B. 

For the full wave calculation, performed in the following, numerical values 
for the atmospheric parameters, like mean pressure p, mean temperature T, 
and velocity of sound C, were taken from CANTOR and COLE (1963). Considering 
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descending waves b * f o at an altitude of z. ~ 31 km (10 mb level), we ask for 
the descending waves b l at z. =» 16 km (100 mb level) under the assumption that 
the ascending waves disappear at this altitude (a 1 « 0). The ascending waves at 
the altitude are then, following Equation (B-3), 

a y = Rj bJ . (15). 

We now specify the descending waves b j either as a gravity wave (bj f 0) 
or a heat conduction wave (bj f 0). From Equations (B-l), (B-4), and (15) the 
physical parameters at the altitude z . follow as 


where 


or 


IV - 


,y^/i 


pij c jV 


G R d 


H^G 

1 

0 


bJ for gravity waves 


R i 
G^H 


h r h j 


b^ for heat conduction waves 


(16) 


(16a) 


(16b) 


Numerical results of Equation (16) are given in Figure 6. Here the following 
values of parameters k x , «, and K eddy from Equations (8), (9), and (13) have been 
used*. 
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k a 1 . 5 x 10” 3 km ” 1 

A 

w = 7 x 10" 6 sec ” 1 

K eddy " 2 x 10 s yW erg cm " 1 sec " 1 deg *' 1 . 


The dispersion equation, Equation (A-2), depends on the Coriolis parameter f . 
However, since the discussed stratospheric warming phenomena occurred only 
over the northern winter hemisphere, and the Coriolis parameter is not of great 
influence on the vertical wave number k z within the considered range of « e f f 
and k x , the parameter f has been chosen constant and the value for 45°N has been 
adopted as representative, in this case. 


Figure 6a gives the amplitude B of the relative temperature disturbance at 
the altitude z., normalized to the value at altitude z. , which is 


(AT/T) 6km 
- 3 e j$ 

(AT/T) 2j = 31kro 


(17) 


as a function of effective angular frequency « eff . Figure 6b shows its phase 0, 
where the right ordinate scale is the vertical travel time of the wave crest which 
is considered to be almost equal to the time delay 


At = - 


0 


CO 


(18) 


The lines marked by n GR" in this figure have been calculated from Equation (16a) 
and are related to an initial gravity wave; the lines marked by n HC n have been 
calculated from Equation (16b) and are related to an initial heat conduction wave. 
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Comparing the data in Figure 5 with the observed values presented in 
Figures 3 and 4, a best fit for an initial heat conduction wave was found for 

co eff = 4.6 x 10' 8 sec" 1 . (19) 

In Figures 3 and 4, the solid lines labelled as I represent theoretical values 
calculated from Equations (16) and (16b). While the theoretical curve in Figure 3 
fits well the observed January 1963 values of time delay of the occurrence of 
temperature maximum with decreasing height, the theoretical (solid) curve I in 
Figure 4 indicates a larger decrease of relative temperature amplitude with de- 
creasing altitude than observed during the 1963 warming. This suggests a larger 
attenuation rate of the temperature disturbance. 

A better fitting pair of curves, the dashed lines labelled as II in Figures 3 
and 4, are theoretical curves calculated with slightly changed wave parameters: 

co = 5 x 10 -6 sec” 1 

k x “ lx 10" 3 km” 1 (20) 

co f f = 2.2 X 10“ 8 sec” 1 . 

e f f 

The assumed value for co is equivalent to an impulse bandwidth of At ~ l/co = 

15 days, which is somewhat larger than the observed value of A r ~ 10 days. 

This discrepancy between observation and theory is not surprising in view of the 
rather primitive theoretical model. In this connection it might be emphasized 
that the theory presented, which is based on a strict perturbation method, neglects 
higher order terms of the perturbation parameters AT or Ap. Thus, the theory 
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cannot predict the heating of surrounding air due to dissipation of the penetrating 
wave. It may, therefore, be expected that the theoretical dissipation rate is 
larger than the real dissipation. These considerations seem to point into the 
right direction in order to explain and correct for the discrepancy mentioned 
above. 

From Figure 5 we see a quasi-linear relationship between time lag At and 
attenuation rate of (AT/T) for heat conduction waves. This relationship can be 
qualitatively verified in Figures 3 and 4. The time lag and the attenuation rate 
of the January 1963 warming were smaller than the corresponding values of the 
February 1952 event, as suggested from Figure 5. 

The numerical value of « ef f in Equation (19) is equivalent to a mean southerly 
wind of 

_ « ~ »ef£ 

u v " : ~ 4.5 x ID” 3 km sec -1 , 

k x 

Since w > > w , , , the exact numerical value of u has not much influence on u . 

err 7 ett x 

Moreover, from Equation (A-2) it follows that « ef f exists only in the product 
« h co f f . Therefore, a change in the numerical value of the austausch coefficient 
k f f , which is inverse by proportional to w h , influences merely the numerical value 
of « ef f but does not significantly change the result of the foregoing discussion. 

In Figure 6, the dashed lines have been calculated from the Wentzel- Kramer- 
Brillouin (WKB-) approximation [Equation (B-6)] . The close agreement of the 
dashed and the solid lines, calculated from full wave analysis in the frequency 
range 10“ 8 < co eff < 1Q“ 7 sec" 1 , in Figure 5 indicates that coupling and reflection 
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are of minor importance in this frequency range and that the nature of the wave 
type does not significantly change during downward propagation. 

In Figures la and lb the dasn-dotted lines are planes of equal phase of the 
downward propagating heat conduction wave calculated from the WKB solution 
[ Kquation (R-8)] using the wave parameters of Equation (20). They agree with 
the observed location of the wave crust of the impulse type disturbance. 

It can therefore be concluded from Figures 1, 3, 4, and 5 that the hypothesis 
of downward propagating heat conduction waves through an atmosphere exhibiting 
a resultant southerly wind in the order of 5 m sec" 1 is consistent with the observed 
wave-like pattern of the three-dimensional temperature field during sudden 
stratospheric warmings in the northern hemisphere in winter. 

6. THE ATMOSPHERE AS A SELECTIVE FILTER WITH RESPECT 

TO ATMOSPHERIC WAVES 

In section 5 it has been shown that descending plane harmonic heat conduction 
waves of angular frequency w ™ 7 x lo~ 6 sec" 1 propagating under a certain angle 
of incidence (horizontal wave number k x = 1.5 x 1CT 3 km" 1 ) are consistent with 
the observed wave like pattern of the temperature field during sudden strato- 
spheric warmings. There remains the question why these waves possess just 
these values of « and k x . 

The impulse shape and the direction of energy flow depend on the impulse 
generator, its antenna characteristics and of the response of the surrounding 
atmosphere to the generated atmospheric wave modes. The determination of the 
impulse generator and its antenna characteristics is outside the scope of this 
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paper. We only can suggest from the observations that the generator creates a 
single impulse with a bandwidth of about 10 days and that the antenna beam is 
mainly directed toward the winter hemisphere. Moreover, we can expect that 
evanescent waves (naving nearly zero real vertical wavenumber) are much less 
excited by the "antenna" than the propagation modes. 

•v 

In order to determine the influence of the atmosphere on the wave propaga- 
tion full wave calculations of the land described in section 5 have been repeated 
for the same atmospheric model (16 £ z £ 31 km) and for downward propagating 
waves. The results are plotted versus horizontal wavenumber k x in Figure 7. 
Here, the parameter is the effective angular frequency « Rff . We notice from 
Figure 7 that the amplitude of the temperature disturbance decreases with in- 
creasing k x and increasing « e f f . The amplitude of the heat conduction wave has 

\ 

dropped to the value £ 10" 2 at w ef f 10" 8 sec" 1 and k x 2 2 X 10' 2 km ' 1 . 

The stratosphere therefore prevents the downward propagation of heat conduction 
waves in these ranges of k x and w e f f . 

With decreasing « e f f and k x the phase angles of the disturbance approach 
the value zero which indicates that in this range heat conduction waves behave 
like evanescent waves. The excitation of evanescent waves by a generator im- 
bedded within the atmosphere is not very likely. Therefore in the ranges 
w ef f < 10* 8 sec' 1 , k x < 10' 4 km* 1 we expect again only a small amount of 
energy from downward propagated heat conduction waves . 

For gravity waves the same arguments are true in an even stronger sense. 
At small horizontal wave numbers (k x < 10" 3 km" *) gravity waves are of an 
evanescent type. At higher values of k x (> 10“ 3 km" 1 ) gravity waves are either 
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attenuated more strongly than heat conduction waves (at w eff < 10" 7 sec' 1 ) or 
their amplitude has dropped by more than 10' 2 (at « ef f > 10" 7 sec" *). 

Therefore we conclude that downward wave propagation within the strato- 
sphere is possible mainly by heat conduction waves with a predominant horizontal 
wave number of the order of k x ~ 10" 3 Ian" 1 and a predominant effective angular 
frequency of w eff ^ 10" 8 - 10' 7 sec" 1 which is consistent with the observations. 
The limitation to a certain value of « eff equivalent to a mean horizontal southerly 
wind of the order of 5 m/sec suggests that it might be this last -condition — 
namely the southerly wind of that strength — which triggers the energy transport. 
If the wind blows from the north, no turbulent heat conduction wave can propagate 
downward into the lower stratosphere. 

The appropriate conditions for the triggering wind system may, perhaps, be 
established by the intensification stage, defined by HIROTA (1967b), of the warm- 
ing phenomenon. 

% - 

7. THE UPWARD TRANSPORT OF WAVE ENERGY 

We finally ask about the upward transport of wave energy. For this reason we 
have calculated in Figure 8 the ratio between the relative temperature disturb- 
ances in z. =70 km and in =45 km height versus horizontal wavenumber k x 
in the same manner as in Figure 7. The variable parameter is again the effective 
angular frequency o ef f . Here the calculations are based on the assumption that 
at the lower height z L a gravity wave (GR) or a heat conduction wave (HC) create 
a temperature disturbance of the amplitude one and that at the altitude z. no 
descending wave exists. 
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From Figure 8 we notice that the amplification of the relative temperature 
disturbance generally does not exceed the factor 1.5. Only gravity waves at 
vertical incidence (k x < 10“ 4 km” 1 ) and low frequencies (co eff < 10" 7 sec" 1 ) 
have higher values. Heat conduction waves with high horizontal wave numbers 
and high frequency ranges (k x > 10" 2 km" 1 ; « ff > 1G~ S sec " 1 ) as well as 
waves with small horizontal wavenumbers and low frequencies (k x < 10“ 3 km" 1 ; 
u ef f < 10" 8 sec" 1 ) can be excluded as participating in an effective upward energy 
transport and there remains a range around u eff ~ 10" 7 sec; k x ^ 10" 3 - 
10" 2 km" 1 for the propagation modes. However, the predominance of gravity waves 
at vertical incidence (k x < 10“ 4 km" 1 ^ in the upward energy transport can not be 
excluded. 

LAUTER (1967) observed that the phenomenon of the winter anomaly is 
increasing with height and that it has perhaps a maximum in the 100 km height 
region. Figure 9 shows the relative temperature disturbance as function of the 
upper boundary z. of our model normalized again by the temperature disturbance 
in 45 km height. The curves marked by n GR” are for gravity waves with the 
parameters o eff = 10" 7 sec" 1 ; k x = 10" 5 km" 1 , those marked by ”HC n are for 
heat conduction waves with the parameters « eff = 10" 7 sec" 1 ; k x = 10“ 3 km’ 1 . 
In these calculations the factor A of the coefficient of heat conductivity has been 
taken as constant and of the value shown in Equation (13), below an altitude of 
85 km. Above 85 km, the value A drops linearly to the value A = 2 x 10 2 erg/cm 
sec deg of the molecular heat conductivity. The amplitudes of both waves show a 
maximum between 80 and 90 km which is due to the decrease of k , , above this 

eddy 

region. The time lag at 65 km height is 2.5 days for heat conduction waves and 4 
days for gravity waves and remains nearly constant below 100 km altitude. We 
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compare this time with the travelling time of downward propagating heat conduc- 
tion waves, considered in Figures 3 and 4, between 45 and 30 km, which is 
about 5 days. 

From correlation studies between sudden stratospheric warmings and winter 
anomalies in the ionospheric D-layei *>t has been found that the maximum of the 
winter anomaly event preceeds the maximum of the stratospheric warming event 
by some days (LAUTER, 1967). During a winter anomaly SECHRIST et al. (1968) 
observed a temperature disturbance within the ionospheric D -region which had 
an impulse shape with maximum amplitude of the order of 50°C. Both observations 

i 

are in general agreement with what we expect from the amplitudes and from the 

} 

fane of phase lag in Figures 7, 8 and 9, though it is, of course, presently 
impossible to draw more than very crude conclusions from a comparison 
between the sparse observational data and the theoretical calculations. 

8. CONCLUSIONS 

The wave-like appearance of the temperature disturbance during the sudden 
stratospheric warming phenomenon of January 1963 can be described by a down- 
ward propagating plane turbulent heat conduction wave of angular frequency 
co = 7 x 10' 6 sec -1 (equivalent to an impulse width of A r = 10 days) and of hori- 
zontal wave length ^ x = 4200 km (equivalent to a horizontal wave number k x = 

1.5 x 10' 3 km" x ). Computations of the horizontal and the vertical phase velocity 
as well as the dissipation rate of the wave amplitude using these wave parameters 
agree with values derived from observation. In these calculations, turbulence 
has been introduced by a turbulent austausch coefficient 10 3 times larger than the 
molecular coefficient of heat conductivity. However, a necessary condition for 
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agreement between theory and observation is the required existence of a mean 
southerly wind of the order of 5 m sec ” 1 , in order to meet the special wave 
characteristics. This mean wind field is suggested to be the triggering mechanism 
of the wave propagation observed during sudden stratospheric warmings,* it might 
be established by processes involved in the intensification period recently re- 
vealed by H1ROTA (1967b). In case of a mean northerly wind, downward propaga- 
tion of waves would be prevented, according to the wave theory. 

From energy and vertical stability considerations it can be concluded that 
the generator for this wave propagation, the production of the initial heat impulse, 
must be located below 50 kilometers, i.e., at or below the stratopause. 

The calculation of upward wave propagation from 45 Ion to 110 km shows 
first a slight increase and then a decrease of the relative wave amplitude, with a 
maximum at 85 Ion. This maximum results from the transition from eddy heat 
transport to molecular heat conduction near the turbopause. The time of arrival 
at an altitude of 75 km is about half the propagation time of downward propagating 
waves between 45 Ion and 30 km the time difference being some days. Both find- 
ings are consistent with observations of the short period winter anomalies in the 
ionospheric D-region. 

Thus, both the observed propagation of sudden stratospheric warmings and 
the observed associated effects of increased D-layer absorption can be understood 
as consequences of a sudden heat impulse released near the stratopause. The 
calculations were based on simple plane harmonic waves; sudden stratospheric 
warmings, however, actually have an impulse structure. Thus, an exact theo- 
retical treatment of the impulse propagation will certainly modify the numerical 
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results. Knowledge of the generation mechanism and of the ,f antenna charac- 
teristics’ 1 of the generator within the stratopause layev will; however, be re- 
quired in order to accomplish this treatment, which therefore, has to be post- 
poned until appropriate observational data become available. 
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Appendix A 


EIGENVALUES OF PLANE TURBULENT ATMOSPHERIC WAVES 

It can be shown (VOLLAND, 1968a) that the coefficient of viscosity ri is 
negligible as compared with the coefficient of heat conductivity k in the treatment 
of the propagation of free internal plane harmonic waves through an isothermal 
atmosphere as long as the conditions 

k 2 c 2 

“h > “a > ~ (A-l) 


hold . It is 


CO 


ef f 


CO 
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= co - k x u x effective angular frequency 
angular frequency 

mean horizontal wind velocity in the direction of wave propa- 
gation (x-direction of a Cartesian coordinate system) 

horizontal wave number 


_ rj= 

w h 2 V 


= LJ 

a 2 C 


y - — ratio between the specific heats at constant pressure and at 

C v 

constant volume 
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g 


gravitational constant 


4 * 


C velocity of sound 

K g 

y = — heat conduction velocity 
c p P 

k coefficient of heat conductivity 

p mean pressure 

T mean temperatu: j. 

Under these circumstances the viscosity waves are decoupled from the 
gravity and from the heat conduction waves. Viscosity waves are heavily attenu- 
ated and do not contribute very much to the energy transport. 

Taking into account the unequalities in Equation (A-l) we find the vertical 
wave number k z of the predominant gravity waves and heat conduction waves 
from Equation (44) in VOLLAND (1968a) 
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with 

f = 2Q sin ? Coriolis parameter 
Q rotational period of the earth 
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geographic latitude 


2/ r - l' 

« ~ — w a Brunt-Vais£U>* frequency 

* I 

The minus sign (v - 1) in Equation (A-2) is due to gravity waves, the plus 
sign (v = 2) is due to heat conduction waves. 

In the treatment of turbulent waves we replace as a first order approximation 
the coefficient of molecular heat conductivity by a coefficient of eddy heat con- 
ductivity: 


K eddy^ LK (L ~ 10 3 ) 

In the range of our interest (10’ 8 < “ e(f < 10 ’ 1 sec" 1 ; z < 100 km) it is 
generally 




w h “off >? J 




k 2 C 2 CO 2 
X g 
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e £ f 


Therefore approximate solutions of Equation (A-2) are 
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(gravity waves) 
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(A-3) 




]/"T co e f£ (1 “ j) (heat conduction waves) 


27 



indicating that in this frequency range turbulent gravity waves are of an evanes 
cent type. For f - o heat conduction waves too become evanescent waves 



Appendix B 


FULL WAVE THEOHY AND RAY APPROXIMATION 


Plane characteristic waves propagating through an inhomogeneous medium 
are coupled with each other. The real atmosphere is inhomogeneous for gravity 
waves and heat conduction waves because the vertical wave number k z is a 
function of altitude via the term 


w h 


yc p p (z) 

K eddy 


Therefore, full wave calculations are necessary to desci*ibe the wave propaga- 
tion through the atmosphere. 


The general solution of full wave analysis has the form (Volland* 1968b) 



(B-l) 


is a column matrix containing the ascending (a 1 ) and descending (b 1 ) gravity 
waves (G) and heat conduction waves (H) in the altitude % and P ij is a 4 x 4 
matrix which connects the characteristic waves c in the altitudes z. and z. . 

i i 
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A rearrangement of the elements of Pi leads to reflection-, transmission-, 
conversion- and coupling coefficients: 
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(B-2) 


and equivalent expressions for R> and P . Now the waves going into a layer are 
connected with the waves coming out of the layer by 




hi 


(B-3) 


They are connected with the physical wave parameters like vertical velocity Aw, 
pressure Ap, temperature AT and temperature gradient d AT/dz by the relation 


where 
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and 


Qiv = (q„ - i A + — ^) f v 
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Hay approximation is valid if the coupling between the different wave modes is 
weak enough so that reflection and conversion as well as coupling can be neglected 
(R ~ 0; h T g - g T h 0). Then the transmission coefficients are approximately 
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With the special normalization of Q [(Equation (B-5)] , Equation (B-6) in connec- 
tion with Equation (B-4), is in fact the WKB-solution. 
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Figure la. Cross-section of Temperature (solid lines) and pressure (dashed lines) along 90°W 
c-jd 90°E on 26 January 1963 derived from stratospheric weather maps (METEOROLOGISCHE 
ABHANDLUNGEN 1963-1966). The dashed-dotted line represents a calculated plane of equal 
phase of heat conduction waves. 
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Figure lb. Cross-section of Temperature (solid lines) and pressure (dashed lines) along 90°W 
and 90°E on 20 January 1963 derived from stratospheric weather maps (METEOROLOGISCHE 
ABHANDLUNGEN 1963-1966). The dashed-dotted line represents o calculated plane of equal 
phase of heat conduction waves. 
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Figure 2. Temperature variciion with time at different isobaric levels at 70 N and 90 W during the January 1963 stratospheric 
warming, derived from stratospheric weather maps (METE0R0L0G1SCHE ABHANDLUNGEN 1963-1966) 
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Figure 3. Time delay A t versus height of the occurrence of the temperature maximum during the 
1963 stratospheric warming (solid circles) derived from Figure 2. Open circles represent cor- 
responding values derived for the February 1952 event, Horizontal bars indicate range of un- 
certainty of the derived data. Solid curve (I) and dashed curve (IX) were computed with the 
following parameters: 

w = 7 x 10“ 6 sec~^ 5 x 10"^ sec~^ 

“ o(f = 4.6 x UT 8 sec -1 2.2 x 10' 8 sec' 1 

k x = 1.5 x lO-Hm' 1 1.0 x TO -3 km -1 
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Figure 4. Relative amplitude of temperature disturbance as function of height. The January 1963 
data (solid circles) were derived from Figure 2. Open circles represent corresponding data for 
the February 1952 event. Horizontal bars indicate range of uncertainty. The solid (I) and dashed 
(II) curves have been calculated with the same parameters listed in Figure 3. The dashed-dotted 
line represents the lower limit of amplitude decrease for a nondissipative wave. 
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Figure 5. Minimum change of temperature amplitude (Curve A), as a function of altitude, required 
to explain a temperature increase of 50°K at 30 km after Equation (12); and maximum permissive 
temperature increase (Curve B) as function of altitude in an isothermal mesosphere, between 
70 km and 50 km, and a colder upper stratosphere (below 50km). The hatched area delineates 
the area of most probable origin of the initial heat impulse. 
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Figure 8. Amplitude and phase of relative temperature disturbance at 70 km altitude as a func 
tion of horizontal wave number k x (in km -1 ) for initial gravity waves (GR) and initial heat con 
duction waves (HC) at 45 km altitude for two different effective angular frequencies co^ (ii 
sec"" 
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